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Abstract 

The axisymmetric deformation and motion of interacting droplets in an im- 
posed temperature gradient is considered using boundary-integral techniques for 
slow viscous motion. Results showing temporal drop motion, deformations and 
separation are presented for equal- viscosity fluids. The focus is on cases when the 
drops are of equal radii or when the smaller drop trails behind the larger drop. For 
equal-size drops, our analysis shows that the motion of a leading drop is retarded 
while the motion of the trailing one is enchanced compared to the undeformable 
case. The distance between the centers of equal-sized deformable drops decreases 
with time. When a small drop follows a large one, two patterns of behavior may 
exist. For moderate or large initial separation the drops separate. However, if 
the initial separation is small there is a transient period in which the separation 
distance initially decreases and only afterwards the drops separate. This behavior 
stems from the multiple time scales that exist in the system. 

Key Words: droplets interaction; thermocapillary migration; drop deforma- 
tion; boundary integral equations 
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1 Introduction 



It has long been known that a viscous droplet submerged into another immiscible non- 
isothermal liquid migrates in the direction of the temperature gradient, VT. This 
droplet's migration is due to the temperature-induced surface tension gradient Vs(7) 
at the droplet interface S and it is addressed as a thermocapillary migration. Young et 
al (|I]) were first to investigate the mechanism of such migration both theoretically and 
experimentally. These authors found conditions for the stabilization of a buoyant bubble 
under gravity force by the thermocapillary force which was directed oppositely. They 
obtained a theoretical prediction for the migration velocity of a single spherical bubble 
of radius a, which is placed in a viscous fluid of viscosity r/i, with an imposed constant 
temperature gradient A. In the limit of zero Marangoni number |^ i}-^-^ with neg- 

ligible convective transport effects) the droplet velocity U in the laboratory coordinate 
frame is related to the uniform temperature gradient by 
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U, U = — -. [1] 
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Here rj = i]2/rii and k = 1^2/1^1 are the ratios of viscosities and thermal conductivities 



dT 



is a constant derivative 



between the internal and the surrounding fluids, respectively, 
of the surface tension with respect to the local temperature T on the interface, and U is 
a dimensionless migration velocity. In the case when the droplet and the external fluid 
have equal viscosity and thermal properties f7 = 2/15 ~ 0.133. A review of the further 
development of the subject can be found in @. 

Bratukhin (|^) considered the migration of a deformable viscous droplet at the small 
Marangoni number. Using perturbation techniques in a spherical coordinate system (r, 6) 
with origin at the center of the droplet, he found an expression for the deformation of 
the droplet's shape r = R{9): 

R{e) = 1 + Ma'^S2P2{0) + 0{Ma^), 

3 (5-2 {l~p){l + v) PrU 49x-7x^ + l8^ \ 
470 V 2 4 + r/ 63x {3 + 2k){2 + k) J ' 



Here Pr = v\lx\ is the Prandtl number, p = P2/P1 and x = X2/X1 are the ratios of density 
and thermal diffusivity between the internal and the surrounding fluids, respectively, 
7o = 1 9x0-1 f]\ is the dimensionless surface tension corresponding to the center of the 
drop, and -P2(^) is the second order Legendre polynomial. 

The flrst term of is related to the inertia term in the Navier Stokes equation, while 
the second one is related to the convective term in the energy transport equation. It is 
interesting to note that the expansion of the droplet deformation begins with the term 
of O(Ma^), and that this term vanishes in the case that combines equal densities of the 
phases (p = 1) and negligible convective heat transport [Pr = 0). It follows from 
that the deformations are small also for large values of the reference surface tension 70 
(small capillary number). 

Balasubramaniam and Chai (|^) showed that the Young solution satisfles the full 
Navier-Stokes equations if the temperature fleld on the spherical interface is proportional 
to cos^, i.e. in the case of negligible convective transport. The perturbation of the 
spherical form constructed in (P under the assumption of small capillary number has 
the same form as 0. 

In real applications, however, fluid drops are not isolated and their motion is influ- 
enced by the presence of neighboring particles. This interaction of drops substantially 
alters the migration velocities and it may cause drops' deformations even in the case of 
negligible convective transport. On the other hand, the second mechanism of the defor- 
mations' suppression (large surface tension) is expected to be valid in the multi body 
case as well. 

Most of the numerous theoretical studies of the interaction of bubbles and droplets 
in the course of their thermocapillary migration were performed under the assumption 
of nondeformable drops (zero capillary number). Meyyappan et al (^ and Feuillbois (P) 
studied the Marangoni induced axisymmetric migration of two bubbles and Meyyappan 
and Subramanian (0) studied the non-axisymmetric case. Acrivos et al (H), Satrape (|) 
and Wang et al ( pUj ) examined ensembles of more than two bubbles. The remarkable 
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finding of tliese studies is tliat a bubble migrates in tlie presence of another bubbles of 
the same size with exactly the same velocity as when isolated. Bubbles of different sizes 
affect the motion of the neighboring bubbles. The influence of the presence of a larger 
bubble on the motion of a smaller one is more pronounced than its influence on the 
motion of the larger bubble. 

These analyses were extended to the case of liquid drops by Anderson (|TTD , Keh and 
Chen |l^, |14D and Wei and Subramanian (|l^). In this case, the migration velocity of 
a drop in the presence of another drops differs from the velocity of an isolated drop even 
when the drops are of equal size. Similar to the case of interacting bubbles, a larger drop 
exerts more pronounced influence on the motion of a smaller one. The thermo capillary 
motion of two non-conducting drops was studied by Loewenberg and Davis (0). 

The literature on the thermocapillary motion of deformable drops is limited in contrast 
to a somewhat related problem concerning the interaction of drops undergoing a buoyancy 
driven motion. For the latter case numerous recent studies revealed a rich variety of 
interaction patterns of deformable drops depending on the Bond number and the initial 



configuration of the system ( see e.g. Manga and Stone (0), Zinchenko et al (0)) and the 
literature cited). Numerical simulations of an axisymmetric buoyancy-driven interaction 
of a leading drop and a smaller trailing drop were recently reported by Davis ([191) . It 
was demonstrated that the trailing drop elongates considerably due to the hydrodynamic 
influence of the leading one. Afterwards, depending on the governing parameters, the 
drops may either separate and return to spherical shape or the trailing drop may be 
captured by the leading one, or one of the drops may break up. 

For the thermocapillary induced motion, the effect of deformability was studied mostly 
by a perturbation technique assuming small deformations. Recently, Rother and Davis 
(^) applied lubrication approximation to study the effect of a slight deformability of the 
interfaces on the thermocapillary-driven migration of two drops at close proximity. 

Zhou and Davis (^O]) presented a recent study in which the drops are free to deform 
in the course of their thermocapillary interaction. They used a boundary-integral tech- 
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nique to study the thermocapillary interaction of a deformable viscous drop with a larger 
traihng drop making no a priori assumptions regarding the magnitude of deformations. 
It was demonstrated that in this case the deformations are small and have a small effect 
on the drops motion at moderate separation. The influence of deformability becomes 
significant only when the drops are close together. Zhou and Davis (pO|) did not adjust 
the surface tension to follow the migration of the drops in the external temperature field 
that would be stationary in the absence of a drop. Hence, the boundary conditions that 
they used were not compatible with the continuous change of the positions of the drops. 

The goal of the present work is to extend the analysis of the thermocapillary interac- 
tion of deformable drops to the cases of equal-sized drops and the leading larger drop and 
to study the influence of the deformations on the relative motion of the droplets. We have 
applied the boundary integral method for the simulations of the thermocapillary induced 
motion of two deformable viscous drops in the case of moderate capillary numbers and 
equal viscosity and thermal properties of the dispersed and continuous phases. 



2 Problem formulation 

Consider two drops embedded in an immiscible viscous fluid with a temperature gradient 
applied along their line of center. In the absence of gravity the drops migrate in the 
direction of the applied temperature gradient one beyond another due to thermocapillary 
forces. The trailing and the leading drop and the continuous fluid are marked by indices 
1, 2 and 3, respectively. All the fluids are assumed to be Newtonian and incompressible. 
The interface tension 7 is assumed to depend linearly on temperature 7(T) = 7(To) — 
|^(T — To). The temperature field approaches a given linear function T = Az + Tq far 
from the dispersed species. The problem is considered in a laboratory coordinate system, 
see Fig. 1. 

Non-dimensional variables are introduced using the following scaling: the radius ai 
of the first droplet for the length, V* = ^ aiA/rj^ for the velocity, ai/V* for the time. 
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A for the pressure and Aai for the temperature. The problem is governed by the 
following dimensionless parameters: the Marangoni number, Ma, the Prandtl number, 
Pr, and the capillary number, Ca, where 



Ma 



^7 



dT 



V3 X3 
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aiA 
7o 



[3] 



The other dimensionless parameters are the ratios of the material properties of the fluids, 
r] = rji/ris, Pi/ps, k = ki/k^ and the geometric parameters R = 02/01 and d = C/ai, 
where C is the distance between the drops along their line of centers. 

We shall be interested in the case of small Marangoni numbers and moderate values 
of the other parameters. Following (|]) one can see that, at the leading order in Ma, the 
inertial and convective transport effects are negligible, and the motion inside and outside 
the droplets is governed by the Stokes equations: 

Vp3 = Av3, X G Q3, 

Vpi = r]A\i, X G fij, i = 1, 2, 

Vv, = 0, xGl]„ « = 1,2,3, [4] 

where Pi and Vj are the pressure and velocity field, respectively. The temperature field 
inside and outside of droplets satisfies the Laplace equation 



ATi = 0, X G fii, i = 1,2,3. 



The continuous fluid is quiescent far from the drops 



V3 ^ 0, as x ^ 00. 



At the interfaces the velocity field is continuous. 



Vo(x)=Vj(x), for X G Ej, i = 1,2, 



and the kinematic boundary conditions require that 



[5] 



[6] 



[7] 



n,: ■ Vi = V„i, for x G S,-, i = 1,2, 
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where Vni is a velocity of the interface Sj in the normal direction. The tangential stress 
balance may be written as 

dT (x) 

rr(n3-ni)-n, = x G S,,, z = l,2, [9] 

and the balance of the normal stresses reads 

n, • (Hs - ni)-ni = - ^.(x)) H„ x G S,, z = 1, 2, [10] 

Here 11^ = — P3I + T3 and Ilj = — pjl + r^Tj with Tj being the rate of deformation tensor, 
Tj is a unit vector tangential to the i*^ interface, and ifj = V ■ rij is the mean curvature 
on the i^^ interface. 

The temperature and the heat flux are continuous 
dT- dT 

Ti = T3, K— ^ = — — for X G Sj, i = 1,2, [11] 
on on 

and at infinity the temperature tends to a given linear function 

T ^ z at |x| ^ 00. [12] 



3 Boundary- integral formulation. Numerical 
method 

The boundary- value problem outlined in the previous section may also be formulated in 
an equivalent integral form (|2lD. Following the classical potential theory ( p2D we can 
obtain a system of boundary-integral equations for the temperature and velocities on the 
interfaces. When the thermal properties and viscosities of the phases are equal for all 
Vti phases the boundary integral equations simplify considerably, and the solution to the 
problem acquires an explicit form. Thus, the temperature and the velocities everywhere 
are given by 

T(x) = z, [13] 
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v(x) = -^/J- 



dS{y). [14] 



Here y G S = Ei + E2 and the kernels for the single layer potentials for Stokes flow are 



with ^ = X — y and ^ = |^|. Note that the dimensionless surface tension equals 1/Ca — z, 
and the natural restriction 7 > implies that unequality z < 1/Ca hold for all points 
on the interfaces during the entire time of the process. Applying this unequality to the 
initial moment t = provides an upper bound for the capillary number 

Ca < ^ -, [15] 

where do is an initial separation between the surfaces of the drops. 



In the axisymmetric case, the integrations over the azimuthal angle in [Q can be per- 
formed analytically, and the problem is reduced to the computation of a one- dimensional 
singular boundary-integral for the interfacial velocity 



V X 



B(x,y)-f(y)c?Z(y). [16] 



Here x = (r^;, Zx) and y = (r, z) are vectors on an azimuthal (r, z) plane in a cylindrical 
coordinate system, a = o"i + cr2 is the union of curves describing the drops' surfaces in 
this plane, and dl{y) denotes the differential arc length of the corresponding contour. 

is a traction jump across the interface. The complete expressions for components of the 



kernel -B(x, y) are given in the Appendix. Note that the kernel |T6| have a logarithmic 
singularity at y ^ x. In the absence of Marangoni effect, when the tangential stresses 
are continuous ( dT/dr = 0), these singularities are often eliminated making use of a 
subtraction technique (^If) , that is based on a well known equality 



/ B(x,y)-n(y)rf/(y) = 0. [17] 



When the tangential Marangoni stresses associated with the temperature gradients on 
the drops surfaces are present, this procedure cannot be simply incorporated into the 
singularity-subtraction formulation. For the purpose of a study of the droplets migration 
it is sufficient to compute the normal component of surface velocity Vn = n(x) ■ v(x). 
We employed a subtraction technique (pT| , p4D in which the normal component of the 



velocity on the interfaces can be evaluated as 

Vn = n(x) ■ v(x) = n(x)- / B(x, y)-f (y)f//(y) 

My) n(x).B(x,y).T(y)rf/(y) 

+ n(x)- / [/„(y)-/„(x*]B(x,y)-n(y)rf/(y), [18] 
No- 
where X* is an arbitrary chosen point in the bulk or on the interface. The choice x* = x 
removes a singularity in the integral. The integrand of the first summand in the right- 
hand side of |T^ tends to a finite limit at y ^ x (see the Appendix) . 

When the separation between the drops surfaces is small, 'near-singularities' appear in 
the integrand if the reference point x is located in the gap region. Following Loewenberg 
and Hinch we facilitate the calculation of integrals at such 'near singular' using 
18| , where x* is the image point across the gap region to x. This approach provides a 



regularization of the 'near-singularity' for the regions on the surfaces that are in close 
proximity. 

In our numerical computations, the drop boundaries in the azimuthal plane are dis- 
cretized using marker points and are approximated by cubic splines. The number of 
points on the two interfaces is equal. The marker points are distributed uniformly with 
respect to the arc length (non-uniform distribution was not needed in our simulation 
because the curvature of the interface remains moderate at all time). The normal and 
tangent unit vectors at the interfaces and the mean curvatures are computed by taking 
the derivatives of the position vector with respect to the arc length. When the normal 
velocities on the interfaces are obtained, the positions of the marker points are advances 
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using Euler or Runge Kutta method. After that the cubic sphne approximation of the 
interfaces is reconstructed and the marker points are redistributed on the new surfaces 
uniformly with respect to the arc length. The Krasny (^) filtering procedure was used 
to smooth the spatial frequency spectrum and to remove noise that was introduced by 
the computation of high derivatives in the code. All Fourier modes below some 'tolerant 
number', which depends on the number of marker points were set to 0. At every step 
the positions of the drops' centers of mass were obtained. The velocities of the centers of 
mass were consequently calculated via numerical differentiation with respect to time. The 
calculations are stopped when the interfaces of the drops intersect or when the surface 
tension at any point on the leading drop becomes negative. 



4 Results and discussion 

We consider the axisymmetric drift of two deformable droplets with different size. For 
simplicity we assume that the droplets and the continuous media have identical viscosity, 
densities and thermal properties. This ensures that the temperature field is continuous 
and has continuous derivatives across the interfaces. The droplets are initially spherical 
and start to move from an initial separation distance. The parameters considered include 
the capillary number Ca, the initial droplets' radii -Ri and R2 and the initial distance 
do between surfaces of the droplets. During the motion we observe the evolution of the 
positions and of the form of the droplets, and the flow patterns inside the droplets and in 
the external fluid. The relative position was measured by the distance between centers 
of mass L{t) and by the distance between the nearest points on the surfaces of the two 
drops d{t). The deformation of a droplet was characterized by the Taylor deformation 
parameter 

D = [19] 

where Az is the length of the drop along the axis of symmetry, and Vmax is the maximum 
radial dimension. Positive and negative values of D correspond to oblate and prolate 

10 



shapes, respectively. 



4.1 Migration velocities. The case of non-deformable drops 

We have tested our numerical code by comparing calculated drift velocities with available 
analytic results. As it was mentioned above, there are analytical results (|ll|, |1^) for the 
thermocapillary induced motion of non-deformable spherical drops in a uniform external 
temperature gradient. These results can be applied for the axisymmetric cases of drift of 
drops with different radii Ri and R2, and for different separation distances between the 
centers of mass. 

For the test of our code we considered a case of small capillary number, Ca = 0.001, 
in which the expected deformation of initial spherical shape of the drops during the 
evolution of the positions is negligibly small. The case of well-separated non-deformable 
drops was consider by Anderson (0). He showed that for droplets of equal radii, the 
velocities are equal and have the form 

U=(l + ^,-l-^)\U\- [20] 



L2 3 

Keh and Chen (13) considered also drops at small distance and derived a general formula 



for the velocities 

Ui = MnlUil + M12IU2I 

U2 = M21IU1I + M22IU2I [21] 

where Mu, M12, M21, M22 are mobility coefficients that were calculated numerically, and 
|Ui|, IU2I are droplets' velocities that each drop would have in the absence of another 
drop. 



The velocities of equal droplets calculated using the formulae [^, |^ and our 
Boundary-Integral code are depicted in Fig. 2 versus the separation distance. The 
good agreement in all distance-limits is visible. In the case of large separations our cal- 
culations, denoted by solid circles, coincide with Anderson's asymptotic formula (curve 
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d), and in the limit of infinite distance between the drops we recover the value of the drift 
velocity 0.133 as in the Young-Bratukhin problem (|T], In the limit of a small separa- 
tion distance our results are in good agreement with the calculations of Keh and Chen 
as well (dashed curve). A comparison for the case of drops of different radii is presented 
in Fig. 3. We observe that the case of the drift velocities of the two non-deformable 
different droplets is in a good agreement with Keh and Chen predictions. 

4.2 Deformation of equal-sized drops. Effects of Ca and veloc- 
ity patterns 

Consider now the case of larger capillary numbers, where considerable deformations are 
expected. The time evolution of the Taylor deformation parameter is shown in Fig. 4 for 
different values of the capillary number and of the initial separation distance. We observe 
two qualitatively different scenarios of the deformation development. For sufficiently 
large initial separations and small Ca (e.g. curves d and e), the deformation is small 
and becomes almost stationary after a short transient period. On the other hand, for 
the relatively small initial separations and large Ca (e.g. curves a and b), the evolution 
of the deformation after the transient period becomes slower, but remains essentially 
non-steady. During the initial period, with time scale of 0(Ca), the deformations occur 
because the non-uniform surface tension is not compatible with the initial spherical form 
of the drop. The motion during this period is driven by the surface tension, rather than by 
its gradient, and the deformation is controlled by the normal component of the interfacial 
stress. At longer time scales the development of the deformations is induced mostly by the 
changes in the mutual positions of the drops (characterized by the separation distance). 
In Fig. 5 the deformation parameter is depicted versus separation distance for equal- 
size drops and different capillary numbers. The initial separation was kept constant at 
Lq = 2.05. The leading drop becomes more deformed than the trailing one. This can be 
explained by the fact that the leading drop is located in a hotter region and, thus, it has 
a lower surface tension. The dependence of well developed deformations on the capillary 
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number is demonstrated in Fig. 6. Here we plotted the deformation parameter of the 
more deformed leading drop in Fig. 5 when the separation distance reached L — 2.04. 
The dependence on small Ca appears to be linear. 

The deviation of the surfaces from the spherical shape \frf(9) + {zi{9) — Z-i)"^ — Ri is 
demonstrated in Fig. 7 for different values of time. One can see considerable deformations 
in the gap region while the outer regions of the drops remain almost spherical. It is evident 
from Figs. 4, 5 and 7 that the leading drop deforms to an oblate shape while the trailing 
one becomes prolate. The mechanism of this deformation is shown in Figure 8 where 
we plot the velocity field for two equal drops at close proximity for time t — 10. The 
reference frame moves with the rear point on the axis of the trailing drop as marked in 
the figure. In this reference frame, the drops are almost at rest except for the gap region. 
The motion is composed of a flow past an aggregate of two drops and the deformation of 
the gap region in the opposite direction. The elongation of the trailing drop to a prolate 
shape and the flattening of the leading drop to the oblate shape are clearly evident. 

Similar to the deformation, the velocities of the drops' centers of mass achieve quasi- 
steady slowly changing values after a short transient period. For two equal drops at 
Ca — 0.2 these are plotted in Fig. 2 versus the distance between the surfaces (points on 
the curves a and c). One can see that, at large separations, the velocities of the two equal- 
sized drops are almost equal and approach the velocities of the non-deformable drops. 
In contrast to this, for sufficiently small separation, the velocities of the two deformed 
drops are different. The leading oblate drop moves slower than the trailing prolate one. 

The time evolution of the surface-to-surface distance d{t) and center-to center distance 
L{t) for Ca — 0.2 and different initial separations is shown in Fig. 9. We observe two 
different behaviors depending on the initial separation. At small initial separations, 
surface to surface distance decreases with time (Figure 9 (e)). For larger separations, 
surface-to surface distance shows a transient behavior in which it grows initially and 
decreases after some period (Fig. 9(b)). This initial period of repulsion increases with 
the initial separation, and for — 0.5 (see Fig 9(c)) the surface-to-surface distance grows 
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during the entire time of the process. To understand this behavior recall that the leading 
drop deforms into an oblate shape, thus, moving apart from the trailing one. The latter, 
in turn, deforms into a prolate shape, however, its deformations are lower due to a higher 
surface tension. As a result, the separation distance increases. The difference in the 
effective surface tension increases with the separation distance, hence, when the distance 
is sufficiently small, this effect is overcome by the simultaneous approach of the droplets' 
centers of mass (Fig. 9 (d) and (f)) and, after a period of repulsion, the drops surfaces 
approach each other. The period of initial 'repulsion' becomes wider for larger initial 
separations, while when the initial distances are small enough, the difference in the rate 
of deformation becomes negligible, the initial repulsion period vanishes (see Fig. 9 e) and 
the drops approach during the entire time of the process. Note that the motion of the 
droplets towards each other is quite slow, and in the cases that we considered drops that 
were initially very close together do not progress to a complete collision. 

4.3 Interaction of Drops of Different Size. 

The interaction of equal size drops is governed by the interplay of two effects: rela- 
tive motion caused by the oblate(prolatc) shape of the leading(trailing) drop, and more 
pronounced deformations of the leading drop in the gap region due to its lower surface 
tension. These effects take place also for unequal drops, but in this case they are combined 
with the relative motion caused by the different sizes. Similar to the case of equal drops, 
after a short transient period the migration velocity of the trailing drop is enchanced 
and the leading drop motion is retarded compared to the case with no deformations. If 
the leading drop is smaller than the trailing one, the separation distance decreases with 
time and the drops come to near contact. For this case, the larger drop remains almost 
spherical, and the smaller one deforms considerably. The form of the drops in close prox- 
imity and the velocity field is shown in Fig. 10 (a) for Ca — 0.2, Ri — 1, R2 — 0.5. 
The reference frame moves with the rear point of the trailing drop. The motion is a 
superposition of a flow past an aggregate of the drops and a motion and deformation of 
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the surface of the leading drop. The larger trailing drop deforms mostly in the near gap 
region. The deformation pattern in this case agrees qualitatively with the one reported 
by Zhou and Davis (pOl). 

The situation is different when the trailing drop is smaller than the leading one. In 
this case the gap between nondeformable drops increases and the drops separate. The 
deformability may result in a more complex interaction between such drops if the initial 
separation is small enough. The velocity field is shown in Fig. 10 (b) for Ca = 0.2, Ri = 
1, i?2 = 1-5. The reference frame moves with the point at the center of the gap region. 
One can observe a very slow motion in the gap and considerable normal velocities at the 
interfaces in the outer regions, indicating the developing deformations. The evolution 
of the separation distances and the deformation parameter are shown in Figs 11, 12 
and 13 for various drop size ratio and Ca. Different curves on each plot correspond to 
different initial separations. It is seen that for the large initial separations, the surface-to- 
surface distance increases during the entire duration of the process, as in the case when 
the deformations are absent. When the initial separation is small the distance is first 
reduced and only increase at the later stages. The transition is controled by the capillary 
number (see Fig. 14) and it seems to occur at separation distance of 0(10^^). 

To understand this unusual behavior, consider the evolution of the deformation factor 
that is shown in these figures. The smaller drop deforms much stronger than the larger 
one assuming a prolate form. This is an anticipated result since the deformation of the 
trailing drop is induced by the pulling presence of the leading one and vice versa. The 
larger drop induces larger deformations. Since in our case the leading drop is larger, 
higher deformations of the trailing drop are expected. At the depicted size ratios this 
effect overcomes the effect of a lower surface tension over the leading drop that has an 
opposite effect in the case of equal drops, where the leading drop is more deformed. 

The deformation of the smaller drop develops very fast during an initial transient 
period and becomes slower afterwards. The interface of a larger drop is initially pushed 
by the trailing one and assumes an oblate shape. These deformations are combined with 
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the relative motion of the drops apart from each other due to the different radii. Unless 
the drops are very close in size, after the initial period the elongation of the trailing 
drop becomes slower, and the distance between the interfaces begins to grow. The drops 
separate and the influence of the trailing drop on the leading one causes its deformation 
to reverse itself and it also becomes prolate (see Figs. 12(c) and 13(c)). When the drops 
are very close in size (see Fig. 11(c)), the evolution their shape is more symmetric and 
resembles the one depicted in Fig. 4 for the equal-size case. 

During the initial period the drops may come very close together and attractive 
London van der Waals forces may results in the rupture of the inter-droplet film and 
cause an eventual coalescence. This process, however, is beyond the scope of our work. 
The distance between the centers of mass of the drops increases monotonically while the 
surface exhibits a transient behavior. During the early stages of the dynamics the gap 
shrinks and the surfaces come to close proximity, while at later stages the drops separate 
as in the case with no deformations. The time in which the initial perturbation of the 
shape is relaxed and the gap acquires a minimum appears to depend linearly on the 
capillary number. In this case the system may result in a higher coalescence rate due to 
diminishing gap in the transient period. 

5 Concluding remarks 

The interaction between two drops in a temperature gradient shows a variety of behav- 
ior patterns. Equal drops always approach each other, however, the surface-to-surface 
distance has a transient behavior during which the separation first grows and later di- 
minishes. When two unequal drops are moving, with the bigger one leading, an opposite 
transient behavior is observed if the initial separation distance is small. Here, the surfaces 
first approach each other and later depart while the drops are moving away from each 
other. This non-trivial behavior takes place during an initial period of the process. Since 
our calculations are based on the creeping flow approximation that neglect inertia of the 
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fluid, these results are applicable only if the flow field in the fluid is established much 
faster than an equilibrium form of the deformed surface. Three different time scales are 
present in the process under consideration: a viscous relaxation time = alp/rj^, capil- 
lary surface relaxation time ts = ?73Cn/7o, and the thermocapillary time t* = rj^/ (^Aj, 




that was chosen as the time scale in our considerations. Since t^/t* = Ma and tg/t* = Ca, 
the condition Ma << Ca is sufficient to ensure that the transient period can be described 
in a creeping flow approximation. 

It is interesting to compare the magnitude of deformation in our case with different 
cases such as buoyancy-driven motion ([T9|) and that resulting from non-linear dependence 
of surface tension on the external field (^). We found that the deformation in the case 
of a thermocapillary migration causes behavior of the drops at close proximity that is 
qualitatively different from that observed in the other cases. The calculated magnitude 
in the thermocapillary was relatively small for both leading and trailing drops and for 
different aspect ratios. The rich variety of interaction pattern typical for a motion driven 
buoyancy is not observed. There are several reasons why it is so. The first is a well-known 
fact that the disturbances of the velocity field, resulting from the thermocapillary-induced 
motion of a single drop, decay with the distance from the particle as This is much 

faster than the 1/r decay in the case of the gravity-driven motion. Thus, the influence 
of the presence of a drop on the motion and deformation of the other one is expected to 
be smaller and to decay more rapidly with the separation distance than the one in the 
gravity-driven motion case. The other important difference of these two processes is that 
the Bond number 



that governs the deformability of the drops undergoing gravity- driven migration, may 
take any values from to infinity, and the interesting behavior is observed when this value 
is relatively large. Bo > 3 (see (|19|)). In contrast to this, the value of the capillary number 
that governs the deformability of the drops in the course of their thermocapillary-driven 
migration is restricted from above if only physically relevant cases of positive interfacial 




Bo 



Apgai^ 



7o 
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tension are considered. In our calculations the no n- vanishing values of Ca were typically 
about 0.2. 

Note that using the linear profile of a surface tension - temperature dependence 
results in restriction on the extent of collective drift of drops. Practically there is some 
moment of time at which the calculation must stop, because after this time the surface 
tension becomes zero at the front edge of the leading drop. In the case of temperature 
dependence this situation implies a phase transition. However, when the surface tension 
depends on concentration of a solute, zero surface tension cannot be achieved. In this 
case the dependence of surface tension on concentration ( for some pairs fiuid-solute) 
can be divided to two different regions. In the first region the dependence of tension on 
concentration is strong and can be described by the linear model considered above. In 
the second region, when the surface tension is diminished, its change with concentration 
is very slow and it can be approximated by a constant (p6D. 

We performed a calculation for the case of equal drops using this simple two-region 
dependence of the interfacial tension letting the linear dependence down the value 
7/ (^Aai^^ = 10~^ and keeping it constant beyond this point. When the drops enter 
the region of low surface tension the deformation pattern changes and more pronounced 
distortion of the spherical shapes are obtained. This deformation pattern is accompanied 
by a clear reduction of the drops' speed and the entire migration is halted. A demonstra- 
tion of such behavior is given in Figure 15. More detailed dynamics of such interactions 
are left for the future study. 
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Appendix: Kernels in Eqn. [16] 

If Tx 7^ 0, i.e. the point x is not on the axis of symmetry 



r \V2 



F{k) + 



m 



Bzr — k 



z — 



m 



rl. — r"^ + iz — z. ^"^ 



z — Zx ( r \ 



k / r \ 



irx + rV + (zx - zY 



Brr = — if) {[rl + r^ + 2{zx-zf]F{k)- 

I X' ^1 X' 



{rx + rf + izx-zf ^^^^1' 
Here F and E are complete elhptic functions of the first and second kind respectively, 

^2 4ra;r 

{rx + rY + {zx - zY' 

and Uz and denote the projection of the normal vector in the azimuthal plane on the 
axis z and r. 

If the point x is on the axis of symmetry, r^; = 0, A; = and the kernels take the form 



B -2Tr ^' + ^(^-^-)' 



, Byz 0, 



R — r) ^^(^ Z^ _ 

Asymptotic expansions as y tends to x 7^ 0): 

= - ln[r^ + - Zxf \ + Czz + ■■■ Brz = Crz + ■ ■ ■ , 

n(x)-B(x, y)-T(y) = c^^n^(x)r^(x) + c^rn;,(x)T^(x) + . . . , 

where Cj-^, c^z, c^^, and Czz are constants. If the curve Uj is parametrized by the arc 
length, the kernels are regular for Vx — ^- 
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Figure legends 



Figure 1. Geometric sketch of two viscous drops immersed in an external temperature gradi- 
ent VT parallel of their axis of symmetry. 

Figure 2. Center of mass velocities for equal drops Ri = R2 = 1- The solid line {d), dashed 
line (6) and the points on these lines are the results of Anderson, Keh and Chen and 
our computations for Ca = 0.001, respectively. The points (a) and (c) correspond 
to leading and trailing deformable drops at Ca = 0.2. 

Figure 3. Droplet velocities for Ca = 0.001 and different aspect ratio. In the case (a, b) 
Ri = 1.0 and R2 = 2.0. In the case (c, d) Ri = 1.0 and R2 = 0.5. The lines b and 
d correspond to the larger drop. 

Figure 4. The evolution of Taylor deformation parameter defined in Eqn. ||T^: (a) Ca = 
0.2, do = 0.1; (6) Ca = 0.2, do = 0.5; (c) Ca = 0.1, do = 0.1; (d) Ca = 0.1, do = 
0.5; (e) Ca = 0.1, do = 2. Positive and negative values correspond to oblate and 
prolate shapes, respectively. 

Figure 5. Deformation parameter for equal-size droplets versus the distance between cen- 
ters of mass: (a) Ca = 0.2, do = 0.05; (6) Ca = 0.1, do = 0.05; (c) 
Ca = 0.05, do = 0.05. Positive and negative values correspond to oblate and 
prolate shapes, respectively. 

Figure 6. Deformation of the leading droplet versus the capillary number at L = 2.04 for 
equal-size droplets. 

Figure 7. The local deformations for trailing (a) and leading (6) equal-size droplets for the 
case Ca = 0.2 and do = 0.01. 1, 2 and 3 correspond to 0.5, 3.3, and 10 units of 
time, respectively. 

Figure 8. Velocity patterns for Ca = 0.2, Ri = 1.0 and R2 = 1.0. The reference frame moves 
with the speed of the marked point. 
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Figure 9. The evolution of the separation and center-to-senter distances for equal-size drops 
at Ca — 0.2 as function of the initial separation (a), {b) Transient dynamics for 
intermediate initial separation, (c, d) Evolution of surface-to-surface separation 
and center-to-center distance for relatively large initial separation, (e, /) Evolution 
of surface-to-surface separation and center-to-center distance for relatively small 
initial separation. 

Figure 10. Velocity patterns for (a) Ca = 0.2, Ri = 1.0 and R2 = 0.5; (6) Ca = 0.15, Ri = 1.0 
and R2 — 1.5. The reference frame moves with the speed of the marked point. 

Figure 11. Interaction between a leading larger drop and a trailing smaller one at close prox- 
imity for Ca — 0.2, Ri — 1.0, R2 — 1.1: (a) evolution of center-to-center distance; 
(b) evolution of surface-to surface separation; (c) evolution of Taylor deformation 
parameter. The numbers 1 to 9 denote the different initial separation. 

Figure 12. Interaction between a leading larger drop and a trailing smaller one at close prox- 
imity for Ca = 0.15, Ri = 1.0, R2 — 1.5: (a) evolution of center-to-center distance; 
evolution of surface-to surface separation; (c) evolution of Taylor deformation pa- 
rameter. The numbers 1 to 6 denote the different initial separation. 

Figure 13. Interaction between a leading larger drop and a traihng smaller one at close prox- 
imity for Ca = 0.12, Ri = 1.0, R2 = 2: (a) evolution of center-to-center distance; 
(6) evolution of surface-to surface separation; (c) evolution of Taylor deformation 
parameter. The numbers 1 to 7 denote the different initial separation. 

Figure 14. The relaxation time tr with minimum separation distance as function of the capil- 
lary number. 

Figure 15. Deformation of equal-size drops migrating into a region of vanishing surface tension. 

The upper plot shows the dependence of surface tension on distance. The lower 
plot demonstrates migration and deformation pattern for various time, (a) t = 0; 

(6) t = 5; {c)t = 19; (d) t = 25. 
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